Successful application of multiscale methods in a real reservoir simulator environment     

K.-A. Lie  O. Møyner  J. R. Natvig  A. Kozlova  K. Bratvedt  S. Watanabe  Z. Li 《Computational Geosciences》2017,21(5-6):981-998

For the past 10 years or so, a number of so-called multiscale methods have been developed as an alternative approach to upscaling and to accelerate reservoir simulation. The key idea of all these methods is to construct a set of prolongation operators that map between unknowns associated with cells in a fine grid holding the petrophysical properties of the geological reservoir model and unknowns on a coarser grid used for dynamic simulation. The prolongation operators are computed numerically by solving localized flow problems, much in the same way as for flow-based upscaling methods, and can be used to construct a reduced coarse-scale system of flow equations that describe the macro-scale displacement driven by global forces. Unlike effective parameters, the multiscale basis functions have subscale resolution, which ensures that fine-scale heterogeneity is correctly accounted for in a systematic manner. Among all multiscale formulations discussed in the literature, the multiscale restriction-smoothed basis (MsRSB) method has proved to be particularly promising. This method has been implemented in a commercially available simulator and has three main advantages. First, the input grid and its coarse partition can have general polyhedral geometry and unstructured topology. Secondly, MsRSB is accurate and robust when used as an approximate solver and converges relatively fast when used as an iterative fine-scale solver. Finally, the method is formulated on top of a cell-centered, conservative, finite-volume method and is applicable to any flow model for which one can isolate a pressure equation. We discuss numerical challenges posed by contemporary geomodels and report a number of validation cases showing that the MsRSB method is an efficient, robust, and versatile method for simulating complex models of real reservoirs.  相似文献   

Implication of Apatite and Anhydrite for Formation of an Iron‐Oxide‐Apatite (IOA) Rare Earth Element Prospect,Benjamin River,Canada          下载免费PDF全文

Mihoko Hoshino  Yasushi Watanabe  Yoshiaki Kon 《Resource Geology》2017,67(4):361-383

The Benjamin River apatite prospect in northern New Brunswick, Canada, is hosted by the Late Silurian Dickie Brook plutonic complex, which is made up of intrusive units represented by monzogranite, diorite and gabbro. The IOA ores, composed mainly of apatite, augite, and magnetite at Benjamin River form pegmatitic pods and lenses in the host igneous rocks, the largest of which is 100 m long and 10–20 m wide in the diorite and gabbro units. In this study, 28 IOA ore and rock samples were collected from the diorite and gabbro units. Mineralogical observations show that the apatite–augite–magnetite ores are variable in the amounts of apatite, augite, and magnetite and are associated with minor amounts of epidote‐group minerals (allanite, REE‐rich epidote and epidte) and trace amounts of albite, titanite, ilmenite, titanomagnetite, pyrite, chlorite, calcite, and quartz. Apatite and augite grains contain small anhydrite inclusions. This suggests that the magma that crystallized apatite and augite had high oxygen fugacity. In back scattered electron (BSE) images, apatite grains in the ores have two zones of different appearance: (i) primary REE‐rich zone; and (ii) porous REE‐poor zone. The porous REE‐poor zones mainly appear in rims and/or inside of the apatite grains, in addition to the presence of apatite grains which totally consist of a porous REE‐poor apatite. This porous REE‐poor apatite is characterized by low REE (<0.84 wt%), Si (<0.28 wt%), and Cl (<0.17 wt%) contents. Epidote‐group minerals mainly occur in grain boundary between the porous REE‐poor apatite and augite. These indicate that REE leached from primary REE‐rich apatite crystallized as allanite and REE‐rich epidote. Magnetite in the ores often occurs as veinlets that cut apatite grains or as anhedral grains that replace a part of augite. These textures suggest that magnetite crystallized in the late stage. Pyrite veins occur in the ores, including a large amount of quartz and calcite veins. Pyrite veins mainly occur with quartz veins in augite. These textures indicate pyrite veins are the latest phase. Apatite–augite–magnetite ore, gabbro–quartz diorite and feldspar dike collected from the Benjamin River prospect contain dirty pure albite (Ab₉₈Or₂–Ab₁₀₀) under the microscope. The feldspar dikes mainly consist of dirty pure albite. Occurrences of the dirty pure albite suggest remarkable albitization (sodic alteration) of original plagioclase (An_25.3–An₆₀ in Pilote et al., 2012) associating with intrusion of monzogranite into gabbro and diorite. SO₄^2? bearing magma crystallized primary REE‐rich apatite, augite and anhydrite reacted with Fe in the sodic fluids, which result in oxidation of Fe²⁺ and release of S^2? into the sodic fluids. REE, Ca and Fe from primary REE‐rich apatite, augite and plagioclase altered by the sodic fluids were released into the fluids. Then Fe³⁺ in the sodic fluids precipitated as Fe oxides and epidote‐group minerals in apatite–augite–magnetite ores. Finally, residual S^2? in sodic fluids crystallized as latest pyrite veins. In conclusion, mineralization in Benjamin River IOA prospect are divided into four stages: (1) oxidized magmatic stage that crystallized apatite, augite and anhydrite; (2) sodic metasomatic stage accompanying alteration of magmatic minerals; (3) oxidized fluid stage (magnetite–epidote group minerals mineralization); and (4) reduced fluid stage (pyrite mineralization).  相似文献   

Nutrients,light and phytoplankton production in the shallow,tropical coastal waters of Bandon Bay,Southern Thailand          下载免费PDF全文

Takashi Yoshikawa  Kosuke Tomizawa  Yuki Okamoto  Kazuya Watanabe  Jintana Salaenoi  Kenichi Hayashizaki  Hisashi Kurokura  Satoshi Ishikawa 《Marine Ecology》2017,38(6)

Phytoplankton primary production and its regulation by light and nutrient availability were investigated in the shallow, tropical coastal waters of Bandon Bay, Southern Thailand. The bay was meso‐eutrophicated and highly turbid, receiving river water discharge. Water column stratification was consistently weak during both rainy and dry seasons. Dissolved inorganic nitrogen (DIN) was higher off the river mouth than in the other regions, suggesting that river water discharge was a main source of DIN. By contrast, dissolved inorganic phosphorus (DIP) showed a significant negative correlation with total water depth, implying that regeneration around the sea floor was an important source of DIP. Surface DIN and DIP showed positive correlations with surface primary production (PP) and water column primary productivity (ΣPP*), respectively. The combined correlation and model analyses indicate that total water depth had an ambivalent influence on water column primary production (ΣPP); shallower water depth induced more active regeneration of nutrients, but it also caused higher turbidity and lower light availability as a result of enhanced resuspension of sediments. Furthermore, there was a vertical constraint for phytoplankton during the rainy season: total water depth tended to be shallower than euphotic zone depth. In conclusion, light limitation and vertical constraint owing to shallow water depth appear to be more important than nutrient limitation for water column primary production in Bandon Bay.  相似文献   

Remote effects of mixed layer development on ocean acidification in the subsurface layers of the North Pacific     

Michio Watanabe  Michio Kawamiya 《Journal of Oceanography》2017,73(6):771-784

Using the outputs of projections under the highest emission scenario of the representative concentration pathways performed by Earth system models (ESMs), we evaluate the ocean acidification rates of subsurface layers of the western North Pacific, where the strongest sink of atmospheric CO₂ is found in the mid-latitudes. The low potential vorticity water mass called the North Pacific Subtropical Mode Water (STMW) shows large dissolved inorganic carbon (DIC) concentration increase, and is advected southwestward, so that, in the sea to the south of Japan, DIC concentration increases and ocean acidification occurs faster than in adjacent regions. In the STMW of the Izu-Ogasawara region, the ocean acidification occurs with a pH decrease of ~0.004 year^?1 , a much higher rate than the previously estimated global average (0.0023 year^?1), so that the pH decreases by 0.3–0.4 during the twenty-first century and the saturation state of calcite (Ω_Ca) decreases from ~4.8 down to ~2.4. We find that the ESMs with a deeper mixed layer in the Kuroshio Extension region show a larger increase in DIC concentration within the Izu-Ogasawara region and within the Ryukyu Islands region. Comparing model results with the mixed layer depth obtained from the Argo dataset, we estimate that DIC concentration at a depth of ~200 m increases by 1.4–1.6 μmol kg^?1 year^?1 in the Izu-Ogasawara region and by 1.1–1.4 μmol kg^?1 year^?1 in the Ryukyu Islands region toward the end of this century.  相似文献   

The effect of a hydrous phase on P‐wave velocity anisotropy within a detachment shear zone in the slow‐spreading oceanic crust: A case study from the Godzilla Megamullion,Philippine Sea          下载免费PDF全文

Katsuyoshi Michibayashi  Tohru Watanabe  Yumiko Harigane  Yasuhiko Ohara 《Island Arc》2016,25(3):209-219

We studied the contributions of plagioclase, clinopyroxene, and amphibole to the P‐wave velocity properties of gabbroic mylonites of the Godzilla Megamullion (site KH07‐02‐D18) in the Parece Vela Rift of the central Parece Vela Basin, Philippine Sea, based on their crystal‐preferred orientations (CPOs), mineral modes, and elastic constants and densities of single crystals. The gabbroic mylonites have been classified into three types based on their microstructures and temperature conditions: HT1, HT2 and medium‐temperature (MT) mylonites. The P‐wave velocity properties of the HT1 mylonite are dominantly influenced by plagioclase CPOs. Secondary amphibole occurred after deformation in the HT1 mylonite, so that its effect on P‐wave velocity anisotropy is minimal due to weak CPOs. Although the HT2 mylonite developed deformation microstructures in the three minerals, the P‐wave velocity properties of the HT2 mylonite are essentially isotropic, resulting from the destructive interference of different P‐wave velocity anisotropy patterns produced by the distinct CPOs of the three constituent minerals (i.e., plagioclase, clinopyroxene, and amphibole). The P‐wave velocity properties of the MT mylonite are influenced mainly by amphibole CPOs, whereas the effect of plagioclase CPOs on P‐wave velocity anisotropy becomes very small with a decrease in the intensity of plagioclase CPOs. As a result, the gabbroic mylonites tend to have weak P‐wave velocity anisotropy in seismic velocity, although their constituent minerals show distinct CPOs. Such weakness in the whole‐rock P‐wave velocity anisotropy could result from the destructive contributions of the different mineral CPOs with respect to the structural framework (foliation and lineation). These results show that amphibole has a high potential for P‐wave velocity anisotropy by aligning both crystallographically and dimensionally during deformation in the hydrous oceanic crust. The results also suggest that the effect of a hydrous phase on P‐wave velocity anisotropy within the detachment shear zone in a slow‐spreading oceanic crust varies depending on the degree of deformation and on the timing of hydrothermal activity.  相似文献   

Enhanced chlorophyll associated with island-induced cyclonic eddies in the eastern channel of the Tsushima Straits     

Goh Onitsuka  Akihiko Morimoto  Tetsutaro Takikawa  Atsushi Watanabe  Masatoshi Moku  Yutaka Yoshikawa  Tetsuo Yanagi 《Estuarine, Coastal and Shelf Science》2009,81(3):401-408

The relationship between island-induced cyclonic eddies and chlorophyll a (chl-a) was investigated using field data and satellite images in the eastern channel of the Tsushima Straits. The maximum chl-a concentration around the leeward eddy of the Tsushima Islands was two or three times greater than that of outside the eddy. Two different mechanisms of chl-a enhancement associated with island-induced cyclonic eddies were found in the post-bloom periods. In summer, when nutrients were depleted in the surface layer, eddy pumping increased the nutrient supply in the euphotic zone, resulting in enhanced chl-a around the shallow thermocline near the eddy core. In late autumn, when the mixed layer deepened over the euphotic zone, the mixed layer depth became shallow due to the doming effect of the cyclonic eddy, therefore, the improved irradiance condition led to an increase in the chl-a concentration in the surface mixed layer. Nighttime satellite visible images showed a number of fishing vessels in the lee region of the Tsushima Islands, implying that the enhanced phytoplankton biomass may have resulted in good feeding conditions for fishes and squids in the Tsushima Straits.  相似文献   

Possibility of recent changes in vertical distribution and size composition of chlorophyll-a in the western North Pacific region     

Hiroshi Ishida  Yutaka W. Watanabe  Joji Ishizaka  Toshiya Nakano  Naoki Nagai  Yuji Watanabe  Akifumi Shimamoto  Nobuhiro Maeda  Michimasa Magi 《Journal of Oceanography》2009,65(2):179-186

Our analysis of the last three decades of retrospective data of vertical distributions and size composition of chlorophyll-a (Chl-a) over the western North Pacific has revealed significant changes of three indices related to Chl-a during summer season, as follows: (1) decreasing linear trend of the proportion of Chl-a in surface layer to that of the whole water column by 0.4 and 2.3% year⁻¹ in the subtropical area along 137°E (STA₁₃₇) during 1972 to 1997 and in the Kuroshio Extension area along 175°E (KEA₁₇₅) during 1990 to 2001; (2) increasing linear trend of the depth of subsurface Chl-a maximum (DCM) by 0.4 and 2.6 m year⁻¹ in STA₁₃₇ and KEA₁₇₅; and (3) decreasing linear trend of larger-size Chl-a (>3 μm) by 0.1 and 2.5% year⁻¹ in STA₁₃₇ and KEA₁₇₅, respectively. Water density (σ _θ ) at 75 m depth had also decreased by 0.006 and 0.05 year⁻¹ in STA₁₃₇ and KEA₁₇₅, respectively. The ratio of biogenic opal to biogenic CaCO₃ in the sinking flux decreased by 0.015 year⁻¹ in the subtropical region from 1997 to 2005. These findings may indicate that the subsurface chlorophyll maximum is deepening and larger phytoplankton such as diatoms has been decreasing during the past decade, associated with the decreasing density of surface water caused by warming in the western North Pacific, especially in the summer.  相似文献   

Seasonal variation of horizontal material transport through the eastern channel of the Tsushima Straits   总被引：2，自引：2，他引：0  

Akihiko Morimoto  Tetsutaro Takikawa  Goh Onitsuka  Atsushi Watanabe  Masatoshi Moku  Tetsuo Yanagi 《Journal of Oceanography》2009,65(1):61-71

We conducted hydrographic observations ten times in the Tsushima Strait to reveal seasonal variations of horizontal material transports such as of heat, freshwater, chlorophyll a, and dissolved inorganic nitrogen (DIN) and phosphorus (DIP) through the eastern channel of the Tsushima Strait (ECTS). The volume, freshwater, and heat transport results are of nearly the same order as results reported in previous studies. The annual mean DIN and DIP transports of 3.59 kmol/s and 0.29 kmol/s are large relative to those of the Changjiang and the Taiwan Strait and are horizontally transported through the ECTS. Nutrient transports are high in July–August and October and low in April and November. Increased nutrient transports in July–August and October are due to the appearance of a cold saline water mass in the bottom layer of the ECTS. Changes in DIN transports in summer and autumn, which account for two-thirds of the total annual DIN transport, would have a large effect on the nitrogen budget and biological productivity in the Tsushima Warm Current region.  相似文献   

Patterns in abundance and size of two deep-water gorgonian octocorals,in relation to depth and substrate features off Nova Scotia     

Shana Watanabe  Anna Metaxas  Jessica Sameoto  Peter Lawton 《Deep Sea Research Part I: Oceanographic Research Papers》2009,56(12):2235-2248

Deep-water corals form unique ecosystems, yet very little is known about factors that regulate their distribution and growth. The abundance and size of two deep-water gorgonian coral species, Paragorgia arborea and Primnoa resedaeformis, and their relationship with depth and substratum cover, were investigated at Northeast Channel, off Nova Scotia, in July 2006. This is the first study to measure abundance and size of these two coral species at depths >500 m in the Canadian Atlantic region. A total of 5 transects between 500 and 1000 m depth were examined using video collected by the remotely operated vehicle ROPOS. Abundance of both species was patchy, but higher at these deeper depths than at <500 m. Abundance generally declined with depth, and was moderately correlated with cover of hard substratum (cobble, boulder, bedrock). These relationships were stronger and less variable for P. resedaeformis than for P. arborea, suggesting that factors such as topographic relief may play an additional role in regulating distributions of P. arborea. Maximum colony height was 125 and 240 cm for P. resedaeformis and P. arborea, respectively, and much greater than recorded for depths <500 m. Overall, colony height and depth relationships were strong for both species, but variable among transects. P. resedaeformis showed a negative relationship with depth, while the opposite was observed for P. arborea, suggesting that the two species are affected differently by factors that vary with depth (e.g. temperature, fishing disturbance). Relationships between colony size and size of attachment stone were stronger for P. arborea, especially for overturned colonies, than for P. resedaeformis, suggesting that availability of suitably coarse substrate may be more important for the long-term persistence of P. arborea colonies.  相似文献   

Estimation of groundwater use by crop production simulated by DSSAT‐wheat and DSSAT‐maize models in the piedmont region of the North China Plain     

Yonghui Yang  Masataka Watanabe  Xiying Zhang  Xiaohua Hao  Jiqun Zhang 《水文研究》2006,20(13):2787-2802

Sustainable water use is in serious crisis in the piedmont region of the Taihang Mountains in the North China Plain, owing to rapid groundwater drawdown. Estimating the water requirement for agriculture, the biggest user of groundwater, will be helpful in understanding groundwater decline. Through the use of DSSAT‐3·5 wheat and maize models, we assessed water use in winter wheat and maize, two staple crops in the region, in 1987–2001. Trends between groundwater change and simulated agricultural water use were compared. The results showed that groundwater decline was sensitive to simulated crop water requirement and irrigation requirement. According to regression analysis, 100 mm of water requirement by cultivated land (mainly wheat and maize) resulted in about 0·64 m of groundwater decline. This relationship might be useful in understanding the regional water balance and to help decision‐makers control groundwater decline through controlling crop water use or through long‐distance water transfer. The study demonstrated the usefulness of using the DSSAT model for estimating crop water use and the effectiveness of clarifying the reason for groundwater decline using the simulation results of water use. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献